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Pion Virtual Compton Scattering (VCS) via the reaction 7r~e — > 7r~e7 was observed in the 
Fermilab E781 SELEX experiment. SELEX used a 600GeV/c ir~ beam incident on target atomic 
electrons, detecting the incident ir~ and the final state tt~ , electron and 7. Theoretical predictions 
based on chiral perturbation theory are incorporated into a Monte Carlo simulation of the experiment 
and are compared to the data. The number of reconstructed events (9) and their distribution with 
respect to the kinematic variables (for the kinematic region studied) are in reasonable accord with 
the predictions. The corresponding n~ VCS experimental cross section is a = 38.8 ± 13 nb, in 
agreement with the theoretical expectation a = 34.7 nb. 
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I. INTRODUCTION 

The electric (a) and magnetic (/3) pion polarizabilities 
characterize the pion's deformation in an electromagnetic 
field, as occurs during "fir Compton scattering. They de- 
pend on the rigidity of the pion's internal structure as a 
composite particle, and are therefore important dynam- 
ical quantities to test the validity of theoretical models. 
Based on QCD chiral dynamics, the chiral perturbation 
theory effective Lagrangian, using data from radiative 
pion beta decay, predicts the pion electric and magnetic 
polarizabilities a, = -(3^ = 2.7 ± 0.4, expressed in units 
of 1CP 43 cm 3 |l] ||. Other theoretical predictions are also 
available |Q. 

The pion polarizabilities are usually investigated via 
their effect on the shape of the measured 77r — ► 'fir Real 
Compton Scattering (RCS) angular distribution, as in 
Ref. jl. Since pion targets are unavailable, pion RCS is 
approximated using different artifices, as shown in Fig. [j]: 
the 7r~ Z — > 7T -Zj Primakoff § and 7p — > "/tt + ti radia- 
tive pion photoproduction reactions || ; or by the crossing 
symmetry 77 — > 7r + 7r~ two-photon reaction H?!. In the 
Primakoff scattering, a high energy pion scatters from a 
(virtual, practically real) photon in the Coulomb field of 
the target nucleus. Values of a measured by these exper- 
iments are given in Table nl They cover a large range of 
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TABLE I. Experimental values of a. 



values and have large uncertainties. New high precision 
pion polarizability measurements are therefore needed. 
Electromagnetic studies with virtual photons have the 
advantage that the energy and three-momentum of the 
virtual photon can be varied independently. In the pion 
7*7r — > 77r VCS reaction, where the initial state photon 
is virtual (far from the quasi-real photons of Primakoff 
scattering) and the final state photon is real, polarizabil- 
ities can be measured in the space like region, inacces- 
sible by RCS [|). We thereby access the so-called elec- 
tric a(q 2 ) and magnetic (3(q 2 ) generalized polarizabilities 
of the pion ||, where ct(0) and /3(0) correspond to the 
RCS a and /3 pion Compton polarizabilities. The q 2 - 
dependent a(q 2 ) determines the change AF(q 2 ) in the 
pion charge form factor F(q 2 ) in the presence of a strong 
electric field. The Fourier transform of a(q 2 ) provides 
a picture of the local induced pion charge polarization 
density [fTpf, Similarly, first experiments and cal- 
culation ]12| have been carried out for proton VCS via 
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FIG. 1. Three pion Compton scattering reactions: (a) 

ir~ Z — > ir~ Zj, (b) 7p — ► 77r + n, (c) e + e~ — > e + e~7v~-K + . 



ep —>■ epj. 

We study experimentally the feasibility of extracting 
the "pion VCS" reaction: 



7T e7, 



(1) 



as a step in developing pion VCS as a new experimental 
tool for pion polarizability measurements. The data were 
taken with the Fermilab E781 SELEX spectrometer jlj]]. 
We used a 600 GeV/c ir~ beam incident on target atomic 
electrons, detecting the incident ir~ and the final state 
7T _ , electron and 7. Theoretical predictions (ljjl^l based 
on chiral perturbation theory are incorporated into a 
Monte Carlo simulation of the experiment and are com- 
pared to the data. With this theory, we calculated the to- 
tal cross section (as described later) of the Eq. ([!]) process 
for a limited kinematic region discussed later [Eq. (0)], 
using the Monte Carlo integration program VEGAS pw . 
The result is a(total) = 34.7 ± 0.1 nb. For the kinematic 
range considered, the integrated cross sections are not 
sensitive to the polarizabilities. We nonetheless chose 
this region in order to obtain sufficient statistics for a 
first study of the reaction. 



II. VCS KINEMATICS AND THEORETICAL 
DIFFERENTIAL CROSS SECTION 

We study reaction |l]) in terms of the following five 
independent invariant variables: 



s = (pi + k) 2 , 

Sl = (k> + q') 2 ,s 2 = (p f +q>) 2 , 
r 2 = ( Pl -p f ) 2 , q 2 = (k'-k) 2 . 



(2) 



Here pi is the 4-momentum of the incoming pion, k is 
the 4-momentum of the target electron, and pf, k' , q' 
are the 4-momenta of the outgoing pion, electron, and 
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FIG. 2. Kinematics of the Eq. (nl) reaction in the labora- 
tory frame. The incoming pion beam with 4-momentum pi 
interacts with the target electron with 4-momentum k and 
produces the outgoing pion with 4-momentum pf, 7 with 
4-momentum q' , and electron with 4-momentum k' . 



P, P f R P 

FIG. 3. The three Feynman diagrams corresponding to the 

7r~e — > 7r~e7. In the one photon exchange approximation, (a) 

and (b) correspond to the BH process, while (c) corresponds 

to the VCS process. Here, q is the 4-momentum of the virtual 

photon, q — k' — k. 



photon, respectively, as shown in Fig. ||. The differential 
cross section of reaction ([!]) in the convention of Rcf. |rijj 
reads as: 



da 



8E f E k >E q/ • fc)2 _ M 2 m 2 (2tt) 5 

\M\ 2 6 i ( Pl + k -p f - k' - q')d 3 Pf d 3 k'd 3 q'. (3) 

The invariant amplitude M contains the complete infor- 
mation on the dynamics of the process. The quantity 
\A4\ 2 indicates the sum over the final states and the av- 
erage over the initial spin states. The fourfold differential 
cross section in terms of independent invariant variables 
of Eq. (^) involves a kinematical function A Jl8j and a Ja- 
cobian matrix A4 defining the phase space of the physical 
areas, and is given by: 



da 



2ml 



ds 1 ds 2 dq 2 dr 2 (2tt) 5 A(s, m 2 e ,ml) 16 (-A 4 ) 1/2 



\M\' 



(4) 



Since the variable s, involving the energies of the beam 
pion and of the target electron, is fixed, differential cross 
section (^) actually depends on four variables. 

In reaction ([!]), the final photon can be emitted ei- 
ther by the electron or by the pion, as shown in Fig. |^. 
The first process is described by the Bethe-Heitler (BH) 
amplitude (Fig. ^|a, b), calculable from quantum electro- 
dynamics. The second process is described by the VCS 
(Fig. |3|c) amplitude. Since the source of the final photon 
emission is indistinguishable, one obtains the following 
form of the matrix element of reaction (|I|) [jl4j : 

\M\ 2 = I ^ MVCS + M BH )(M VCS * + M BH *) = 



\M 



BH\2 



\M 



VCS\2 



\ m vcs m bh* + m vcs *M bh \, (5) 



where A4 BH and A4 VCS are amplitudes of the BH and 
VCS processes. 

The general features of the four-fold differential cross 
section can be inferred from Eq. (Q) and matrix element 
calculations. The si-dependence is dominated by the 
(si — ml)^ 1 pole of BH, the cross section varies accord- 
ingly, and is only slightly modified by the s% dependence 
of the VCS amplitude. The S2 dependence is dominated 
by the (s2 — M 2 )^ 1 pole of the VCS amplitude with mod- 
ification of S2 dependence of the BH amplitude, but in 
this case the modification is not as small as in the case 
of si. The r 2 dependence is determined by the r pole 
of the BH amplitude, and the q 2 dependence follows the 
q~ 4 pole behavior of typical electron scattering. The en- 
ergy of the outgoing pion is expected to be high while 
the angle is expected to be small according to the r~ 4 
behavior of the cross section. The energy of the outgo- 
ing photon is mainly expected to be low, as follows from 
the infrared divergence of the BH amplitude. The an- 
gular behavior of the outgoing photon is determined by 
the (sx — to 2 ,) -1 and (s 2 — Af 2 ) -1 poles of the BH and 
VCS amplitudes. The more interesting photons related 
to the generalized polarizabilities are expected to have 
higher energies. The behavior of the outgoing electron is 
completely determined by the g -4 behavior of the cross 
section. 



III. EXPERIMENTAL APPARATUS AND 
TRIGGER 

Our data were taken with the Hadron-Electron (HE) 
scattering trigger Q of experiment E781/SELEX at 
Fermilab. SELEX uses a negatively charged beam of 
600 GeV/c with full width momentum spread of dp/p=± 
8%, and an opening solid angle of 0.5 /isr. The beam 
consisted of approximately 50% it~ and 50% SE- 



LEX used Copper and Carbon targets, totaling 4.2% of 
an interaction length, with target electron thicknesses of 
0.676 barn -1 and 0.645 barn -1 , respectively. 

The experiment focused on charm baryon hadropro- 
duction and spectroscopy at large-x^- The spectrometer 
hosted several projects which exploit physics with a small 
number of tracks compared to charm. SELEX had good 
efficiency for detecting all particles in the final state since 
the produced particles and decay fragments at large- a; f 
production are focused in a forward cone in the labora- 
tory frame. Other requirements in the charm program 
for background suppression include good vertex resolu- 
tion and particle identification over a large momentum 
range. 

Four dipole magnets divide SELEX into independent 
spectrometers (Beam, Ml, M2 and M3) each dedicated 
to one special momentum region. Each spectrometer had 
a combination of tracking detectors. The Ml, M2 and 
M3 spectrometers included electromagnetic calorimeters. 
The 7r — e separation for hadron-electron scattering used 
the M2 particle identification transition radiation detec- 
tor and also the electromagnetic calorimeter. 

The HE scattering trigger was specialized for separa- 
tion of hadron-electron elastic scattering events. The 
trigger used information from a charged particle detec- 
tors just downstream of the target and, from a hodoscope 
just downstream of M2, to determine charged particle 
multiplicity and charge polarity. For the HE require- 
ment, no electromagnetic calorimeter information was 
included. Therefore, the data collected via this trigger 
include hadron-electron elastic and inelastic scattering 
events. 

IV. MONTE CARLO SIMULATION 

Monte Carlo simulations were carried out for tt~ VCS 
signal (0) and background event distributions with re- 
spect to the four invariants Si, Sa, q 2 , r 2 . We used the 
SELEX GEANT package GE781 @. The Monte Carlo 
(MC) study was carried out in four steps: 

1. A VCS event generator was written to search for 
the regions of phase space where the data are sen- 
sitive to pion structure. Several event generators 
were made to simulate a variety of expected back- 
grounds. 

2. The event generators were implemented into the 
simulation package. We studied the resolution, 
detection efficiency, geometric and trigger accep- 
tances for the signal and background. 

3. The offline analysis procedure was developed and 
tuned to devise software cuts eliminating back- 
ground while preserving a VCS signal. 



4. Finally we estimated the expected number of tt 
VCS events. 

The VCS event generator is written, based on differen- 
tial cross section (|]), matrix element calculations, and 
3-body final state kinematics. The acceptance-rejection 
method |^0| is used for event generation. The VCS cross 
section increases rapidly when the direction of the outgo- 
ing real photon is close to the direction of one of the out- 
going charged particles (due to the (si— f"- 2 ) -1 pole of BH 
and the (s2 — M 2 ) -1 P°l e of VCS), or when the energy of 
the outgoing real photon comes close to zero (due to the 
infrared divergence of the BH ((si— m 2 — q 2 +r 2 )~ 1 pole). 
Therefore, if events are generated in the pole region, then 
the efficiency of the acceptance-rejection method can be 
rather low, for the more interesting non-pole regions. In 
order to generate events at an acceptable rate, the pole 
regions are eliminated. Invariants are generated in the 
following regions: 

1000m 2 < si < M 2 , 
2M 2 <s 2 < M 2 , 
-0.2 GeV 2 < q- < -0.032 GeV 2 , 
-0.2 GeV 2 < r 2 < -2m e E^(min) + q 2 + s 1 - m 2 . (6) 

For the photon minimum energy, we choose Ej (min) = 5 
GeV to cut the infrared divergence of BH, and to be 
above the calorimeter noise. Since the VCS calculation 
does not explicitly include the p resonance, we choose up- 
per limits of M 2 for s\ and In Fig. [|, we show the gen- 
erated distribution of events plotted with respect to the 
Mandelstam invariants, without correction for any accep- 
tances. The solid and dashed curves are for BH+VCS 
and BH respectively. The VCS amplitude clearly affects 
the shape of these distributions for Si and S2- For q 2 and 
r 2 the effect is difficult to see for the statistics shown. 
Taken together, the experiment therefore is potentially 
sensitive to the pion VCS amplitude. Since 7r/e separa- 
tion is not 100% efficient, all interactions which produce 
two negatively charged tracks and a photon in the final 
state can generate the pattern of the pion VCS; i.e., can 
create background to the required measurement. For the 
background simulation, as well as in case of the VCS 
simulation, we require 5 GeV minimum energy for the 
photon. 

The dominant background process for pion VCS is ir~e 
elastic scattering followed by final state interactions of 
the outgoing charged particles, such as Bremsstrahlung. 
The MC simulation shows that 28% of the original 
7r~e elastic scattering events generate more than 5 GeV 
Bremsstrahlung photons somewhere in the SELEX ap- 
paratus. We therefore consider the S3 invariant mass of 
the outgoing tt~ e system: 

S3 = ( Pf + k') 2 - (7) 
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FIG. 4. Generated distributions of events plotted with re- 
spect to the Mandelstam invariants(solid line corresponds to 
the generation with full VCS amplitude and dashed line corre- 
sponds to the generation where only BH amplitude was used. 



We expect S3 — s for elastic scattering, and S3 = 
si — s 2 + ml + M% for VCS. Fig. || shows generated and 
reconstructed elastic and VCS events subjected to the 
same charged particle track reconstruction and trigger 
requirements as the data. The data shown are only in 
the kinematic region of Eq. As seen from Fig. || the 
events in the data at low S3 do not arise from ir~e elas- 
tic or 7T _ VCS. To describe the events with low S3, we 
simulate backgrounds (with corrections for acceptances) 
from the interactions: 



M e, 



tt" A(Z) -> M A(Z), 



(8) 



(9) 



where M in Eqs. (g) and (^|) is an intermediate meson 
state which can decay via 7r~7r°, 7r~ 77, 7r~w, etc. Consid- 
ering the threshold energies of reaction (||) , only p meson 
production is allowed at SELEX energies. Fig. |3| shows 
the simulated S3 distributions for reactions (pfjfiP, with 
the same acceptance requirements as in Fig. |T 

SELEX GE781 allows us to estimate the background 
rate from ir~ e elastic scattering. However, it is very diffi- 
cult to estimate the background rate from all neutral me- 
son (tt , etc.) production reactions, such as those shown 
in Fig. ||. The data of Fig. ^(bottom) are qualitatively 
well described by a combination of elastic (Fig. || top) 
and meson production (Fig. ^ bottom) channels. The 
contributions of VCS (Fig. || middle) and p production 
(Fig. ^| top) are relatively much lower. We do not show 
the quantitative sum of all contributions, because of the 
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FIG. 5. Comparison of simulated S3 distributions for 7r~e 
elastic scattering (top) and pion VCS (middle) with the data 
(bottom). 



difficulty to estimate the absolute yields of all meson pro- 
duction channels. Instead, we seek a set of cuts which 
remove as "completely" as possible the background from 
all 7 sources, arising from neutral meson decays. 

Since we measure all outgoing particles, the reaction 
kinematics are overdetermined. Therefore, in the data 
analysis, a constrained \ 2 fitting procedure |2|] is used. 
A veto condition is used based on a 2-body final state 
constrained % 2 kinematic fit for reduction of the back- 
ground from 7r~e elastic scattering p2[ . A 3-body final 
state constrained x 2 kinematic fit is used for extracting 
the pion VCS signal. 

A final state electron can arise from photon conver- 
sion. For this background, the electron is not created 
at the same vertex as the pion. Consequently, the qual- 
ity of the 3-track vertex reconstruction should be low. 
Also, the simulation shows that no VCS outgoing parti- 
cles hit the upstream photon calorimeter. Therefore, in 
addition to the constrained kinematic fitting cuts, we use 
restrictions on the vertex quality and on the total energy 
deposit in the upstream photon calorimeter to suppress 
backgrounds. 

An additional way to reduce the background from n~e 
elastic scattering is to use the invariant mass si and 
the Q Sl angle (angle between pi and (q' + k'); see Fig. ||). 
The simulation shows that if the final photon is emit- 
ted via electron bremsstrahlung, the value of si should 
be low. On the other hand if the photon and electron 
are produced via 7r°, 77, or ui meson, si will "remem- 
ber" the mass of the parent particle. Holding the value 
of the si invariant to be between the squared mass of 
7T° and 77 mesons cuts the background from the electron 
bremsstrahlung and from reactions (||). Fig. [?] shows the 
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FIG. 6. Background via meson production: simu- 
lated S3 distribution for n~e — > pe' (upper), and for 
7r~ A(Z) — > MA(Z) (lower), with arbitrary normalization. 
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FIG. 7. sfsx, ©a! correlation for (a) VCS, (b) 7r~e elastic 
scattering, (c) intermediate meson, and (d) p meson produc- 
tion. Only events inside the O sl -y / sT region indicated by the 
box were accepted for the further analysis. 



distributions of VCS and background simulations. It is 
seen that the region with higher y/si and lower Sl are 
mostly populated with pion VCS. Another possibility to 
reduce the background from interactions (|[^) is to cut 
on the invariant M, defined as: 

M = y/{P^+P e -P e ,y. (10) 

For elastic events M = M w , for VCS M — s/sE, for p 
production M = M p (see Fig. ||). The set of the final 
cuts follows: 

• 1, fulfillment of the pion VCS pattern with: iden- 
tified electron; E 7 >5 GeV for photons observed at 
laboratory angles less than 20 mrad in the down- 
stream electromagnetic calorimeters; no additional 
tracks and vertices; the total energy deposit is 
less than 1 GeV in the upstream electromagnetic 
calorimeter, covering detection angles greater than 
30 mrad. 

• 2, Eq. (JsJ) ranges for the invariants, with 
si(min)=0.0225 GeV 2 . 

• 3, xtiastic > 20, following a constrained fit proce- 
dure g@. 

• 4, Xycs — 5> following a constrained fit proce- 
dure 00. 

• 5, 6 Sl <2 mrad, M < 0.625 GeV. 

To estimate the number of expected VCS events, as well 
as the yields of other n~e elastic or inelastic scattering 
events, we use: 



Nne = N n ■ G ■ Nt ■ (-ex ■ £r ■ (11) 

Here N ve is the number of events observed for a particu- 
lar 7re reaction, N n is the number of incident beam pions 
(~ 4.4 • 10 10 as measured by beam scalers and includ- 
ing particle identification), a is the cross section for the 
particular reaction, e r is the offline reconstruction effi- 
ciency (36.6% for elastic, 2.65% for VCS) of the studied 
process, and Nt is the target electron density. Since not 
all experimental properties are implemented in GEANT, 
an additional efficiency factor e ex is included in Eq. 
This efficiency factor is common to 7re elastic and pion 
VCS reactions. The value e e x is calculated by compari- 
son of the actual number of observed 7re elastic scattering 
events to the expected number of events. The common 
efficiency arises since these two reactions have practically 
the same q 2 dependence; q 2 being the only kinematical 
parameter relevant for the trigger and first order data 
reduction procedure. We calculate the experimental effi- 
ciency e ex (13.4%) from ire elastic scattering, as described 
in Ref. [|14||. For extraction of the reference ne elastic 
scattering events, we employ the cuts of the SELEX ire 
elastic scattering analysis |22j| . The cuts described above, 
considering the studied sources of background, improve 
the signal/noise ratio from less than 1/400 to more than 
361/1. For these estimates, we used the following cross 
sections: ayes = 34.7nbarn and <r e ; astic = 4.27/ibarn 
for 7re scattering; <Jp r i m akof f(C target) = 0.025 mbarn 
and <Jprimakof f(Cu target) = 0.83 mbarn for the in- 
dependent Primakoff scattering. Based on the calculated 
pion VCS cross section, we expect ~8 events in the 7r - 
data sample. 



6 









(a) 


5f...i..ti,,, 


L. 




0.2 0.4 0.6 0.8 1 

M [GeV] 









(t>) 


T T I T T , I T , , 1 





M [GeV] 



1400 r 

1200 ~ r 

1000 r 

800 r 

600 - 

400 r 

200 r 

Ej-l 



(d) 



0.6 0.8 1 

M [GeV] 



FIG. 8. Generated distribution of the M invariant for: (a) 
VCS, (b) elastic scattering, (c) meson production reactions 
([]), and (d) p production reaction (^). The vertical lines 
show the 0.625 GeV M-cut position. 



We consider data only in the kinematic region of Eq. [|. 
On the accepted data set, we applied the system of the 
cuts discussed above. The working statistics with these 
cuts on the data are given in Table III and in Fig ^, where 
we show the effects of the cuts on the S3 invariant. The 
effect of the cuts on the data (Table p^ ) is comparable 
to the effect on the simulated VCS events (Tabic [n]). 



cuts 



% of remain events 



26.8 
13.9 
9.93 
0.61 
0.13 



TABLE III. Percentages of events remaining after using 
the cuts. 



The effect of the above enumerated cuts on the VCS 
signal, and on the backgrounds coming from 7re elastic 
scattering and Primakoff meson production, are listed in 
Table |l[ The results are based on the estimated relative 
cross section for these three processes H]. The cuts are 
very effective in removing backgrounds while retaining 
signal events, such that the final event sample is essen- 
tially pure pion VCS. 



Cuts 


VCS 


Elastic 


Meson production 


1 


29.8 


32.1 


0.24 


2 


9.97 


9.21 


0.03 


3 


8.89 


0.03 


0.03 


4 


3.58 


0.003 


0.004 


5 


2.56 


< 0.001 


< 3. x 10~ 6 



TABLE II. Percentages of the remaining events after using 
the cuts for the MC simulated tv VCS and background events. 




S 3 [GeV 2 ] S 3 [GeV 2 ] 




S 3 [GeV 2 ] S 3 [GeV 2 ] 
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S3 [GeV 2 ] S, [GeV 2 ] 

FIG. 9. Effect of the cuts on S3 invariant, (a) raw distribu- 
tion, (b) after cut 1, (c) after cut 2, (d) after cut 3, (e) after 
cut 4, (f) final distribution of S3 invariant. 



V. DATA ANALYSIS 

In the first stage analysis, events containing one iden- 
tified electron are selected. For these events, particle 
trajectories are checked if they form a vertex inside the 
target material. An event is accepted if it contains ex- 
actly three tracks, including the beam particle and an 
electron candidate, and forms one vertex in the target. 



Finally 9 events (with a statistical uncertainty ±3) were 
accepted as pion VCS. The corresponding 7r VCS experi- 
mental cross section based on Eq. ( |Tl| ) under the assump- 
tion that the background has been completely eliminated 
is a = 38.8 ± 13 nb, in agreement with the theoretical ex- 
pectation a — 34.7 nb. The error given is only statistical, 
and does not include possible systematic uncertainties in 
the efficiency product e ex ■ e r in Eq. (pd|). 

The comparisons of reconstructed (data) and gener- 
ated (theory) event distributions, normalized to unit area 
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are shown in Fig. |10| with respect to the four invari- 
ants Si, S2, q 2 and r 2 of Eq. (0), shown with binning 
that matches the experimental resolutions. The resolu- 
tion for each variable was found by Monte Carlo simula- 
tion, comparing reconstructed and generated events. To 
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MC 
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0.125 
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0.025 



(b) 



. data 
MC 



s, [GeV 2 ] 
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data 
MC 



0.9 
0.8 


L Id) 




0.7 


i data 




0.6 


L... MC 




0.5 






0.4 






0.3 






0.2 






0.1 





-0.2 -0.15 -0.1 -0.05 

q 2 [GeV 2 ] 



-0.2 -0.15-0.1 -0.05 

r 2 [GeV 2 ] 



FIG. 10. Comparison of data and theory normalized distri- 
butions with respect to: (a) si, (b) S2, (c) q 2 , and (d) r 2 . The 
solid line corresponds to data, the dashed line corresponds to 
simulations with a = 2.7. The simulation with a = 6.8 gives 
practically the same result. 



check whether the data and theory (MC) distributions 
are consistent with each other, we use the K-S test of 
Kolmogorov and Smirnov |25|. 



The K-S test is based 
on normalized cumulative distribution functions (CDF). 
We use the K-S D statistic, as a measure of the overall 
difference between the two CDFs. It is defined as the 
maximum value of the absolute difference between two 
CDFs. The significance level for a value D (as a dis- 
proof of the null hypothesis that the distributions are 
the same) is given by the probability P{D) p5fl . A high 
value of P{D) means that the data and theory CDF are 
consistent with one another. 

Following the K-S procedure, we calculate the normal- 
ized cumulative distributions of data and theory, corre- 
sponding to Figs. [To|a-d.The K-S D statistic, and the K-S 
probabilities for consistency of theory/data distributions, 
are given in Table and Figs. 0- |l|. The experimen- 
tal and theoretical CDFs for Si and r 2 look similar. For 
S2 and q 2 , some regions of q 2 (S2) have the experimental 
CDFs larger (smaller) than the theoretical CDFs. The 
values of P(D) (see Table IV) are sufficiently high for 
all five CDF's, as expected if the experimental data are 
consistent with theory. In addition, the prediction of a 



variable 


D 


P(D) 


Sl 


0.18 


0.90 


S2 


0.27 


0.83 


S3 


0.27 


0.92 


q 2 


0.18 


0.99 


r 2 


0.07 


0.99 



TABLE IV. Value of K-S D statistic, and probabilities P, 
for comparison of data with theory for a = 2.7. The compari- 
son of data with theory for a = 6.8 gives practically the same 
result. 



total of 8 events is in agreement with the observed 9 ± 3 
data events. This further supports the conclusion from 
the K-S test that we observe pion VCS events. 

From the limited statistics and sensitivity of this first 
pion VCS experiment, we cannot determine the a polar- 
izability value, nor can we determine which value of a is 
preferred. In a future experiment, the sensitivity to pion 
polarizability may be increased jL4| by achieving a data 
set in which the final 7 (tt) has higher (lower) energies. 
However, such data correspond to a lower cross section, 
and therefore require a high luminosity experiment. 



VI. CONCLUSIONS 

The pion Virtual Compton Scattering via the reaction 
7re — > n'e'j is observed. We developed and implemented 
a simulation with a VCS event generator. We defined 
cuts that allow background reduction and VCS signal 
extraction. The measured number of reconstructed pion 
VCS events, and their distributions with respect to the 
Mandelstam invariants, are in reasonable agreement with 
theoretical expectations. The corresponding tt VCS ex- 
perimental cross section is a — 38.8±13nb, in agreement 
with the theoretical expectation a — 34.7 nb. 
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gives practically the same result. 
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